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We report on a search for direct scalar bottom quark (sbottom) pair production in pp collisions 
at v^ = 1.96 TeV, in events with large missing transverse energy and two jets of hadrons in the 
final state, where at least one of the jets is required to be identified as originating from a b quark. 
The study uses a CDF Run II data sample corresponding to 2.65 fb^^ of integrated luminosity. The 
data are in agreement with the standard model. In an R-parity conserving minimal supersymmetric 
scenario, and assuming that the sbottom decays exclusively into a bottom quark and a neutralino, 
95% confidence-level upper limits on the sbottom pair production cross section of 0.1 pb are obtained. 
For neutralino masses below 70 GeV/c^, sbottom masses up to 230 GeV/c^ are excluded at 95% 
confidence level. 
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Supersymmetry (SUSY) [1[ is an extension of the stan- 
dard model (SM) that naturally solves the hierarchy 
problem [2| and provides a possible candidate for dark 
matter in the Universe. SUSY doubles the SM spectrum 
of particles by introducing a new supersymmetric part- 
ner (sparticle) for each particle in the SM. In particular, 
a new scalar field is associated with each left- and right- 
handed quark state, and two SUSY mass eigenstates qi 
and 92 result from the mixing of the scalar fields. In some 
SUSY scenarios, for which the ratio of the vacuum expec- 
tation values of the two supersymmetric Higgs fields be- 
comes large, a significant mass difference between eigen- 
states in the sbottom^uark sector can occur, leading to a 
rather light sbottom 61 mass state. In a generic minimal 
supersymmetric extension of the SM (MSSM) and as- 
suming R-parity conservation [l| , sparticles are produced 
in pairs and the lightest supersymmetric particle (LSP) 
is stable and identified as the neutralino Xi- Assuming 
a SUSY particle mass hierarchy such that the sbottom 
decays exclusively as 61 — >■ bxi, the expected signal for 
direct sbottom pair production is characterized by the 
presence of two energetic jets from the hadronization of 
the bottom quarks and large missing transverse energy 
$!rp [3| from the two LSPs in the final state. Results on 
searches in this channel using Tevatron data have been 
previously reported by both the CDF and DO experi- 
ments in Run II [J| [5|. This Letter presents new results 
based on an almost ten times larger data sample collected 
by the CDF experiment, corresponding to 2.65 fb~^ of 
integrated luminosity. 

The CDF II detector is described in detail else- 
where ^]. The detector has a charged particle tracking 
system that is immersed in a 1.4 T solenoidal magnetic 
field coaxial with the beam line, and provides coverage 
in the pseudorapidity [3] range \r]\ < 2. Segmented sam- 
pling calorimeters, arranged in a projective tower geome- 
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try, surround the tracking system and measure the energy 
of interacting particles for \ri\ < 3.6. Cherenkov counters 
in the region 3.7 < jryj < 4.7 measure the number of 
inelastic pp collisions to determine the luminosity [7[ . 

The data are collected using a three-level trigger sys- 
tem that selects events with at least one calorimeter clus- 
ter with Et above 100 GeV. The events are then re- 
quired to have a primary vertex with a z position within 
60 cm of the nominal interaction point. Jets are re- 
constructed from the energy deposits in the calorimeter 
towers using a cone-based algorithm [81 with cone ra- 
dius R = yA^ 



A2 - 



0.4, and the measured E^j 



is corrected for detector effects and contributions from 
multiple pp interactions per crossing at high instanta- 
neous luminosity, as discussed in Ref . [9| . The events are 
initially selected with $rp > 10 GeV and two jets with 
corrected transverse energy Ei^ > 25 GeV and pseudo- 
rapidity |7yJ°'| < 2.0, where at least one of the two jets is 
required to have \if'^*^\ < 1.1. Events with additional jets 
with e'^^ > 15 GeV and |7yJ<=*| < 2.8 are rejected. The 
effect of the trigger selection is studied using control data 
samples, parametrized as a function of $!j, and the trans- 
verse energy of the leading jet E'r^ , and consistently 
taken into account [10l|. At a given $!rp, the trigger effi- 
ciency increases with increasing _E;^ ^ ' . Finally, at least 
one of the two leading jets is required to originate from 
a 6-quark jet candidate, using the default CDF 6-jet tag- 
ging algorithm (secvtx) [11|, based on the presence of 
a displaced vertex due to the decay of a 6 hadron inside 
the jet. 

The SM QCD multijet background contribution with 
large $Jj,, due to the mismeasurement of the jet ener- 
gies in the calorimeters, is suppressed by requiring an 
azimuthal separation A(j){]^rp —jet) > 0.4 for the two jets 
in the event. The SM background contributions with 
energetic electrons [ij] from W and Z decays and re- 
constructed as jets in the final state are suppressed by 

■jet 
T,cm 



requiring £'j?om/^T'^ < ^-^ fo'^ ^^"^^ J^t> where E; 
denotes the electromagnetic component of the jet trans- 
verse energy. In addition, events with isolated tracks 
with transverse momentum pr.track above 10 GeV/c are 
vetoed, thus rejecting backgrounds with W ox Z bosons 
decaying into muons or tau leptons. Beam-related back- 
grounds and cosmic rays are eliminated by requiring an 
average jet electromagnetic fraction /om > 0.15 and an 
average charged particle fraction /ch > 0.15, as defined 



13| . The requirements on /em and /ch reject events 



with anomalous energy deposition in the hadronic section 
of the calorimeter or energy deposits in the calorimeter 
inconsistent with the observed activity in the tracking 
system, and have no significant effect on SUSY signal 
and SM background physics samples. 

The dominant source of background in the analysis is 
due to events with a light-flavor jet which is misiden- 
tified as a 5-jet (mistags). There are also contributions 



from heavy-flavor QCD multijet events with large Ifjj. and 
passmg the A(/)(_^y— jet) requirement discussed above. In 
both cases, the background estimation is extracted from 
data using a procedure similar to that described in [ll] 
and p^], where different data-driven multidimensional 
parameterizations are employed to estimate the proba- 
bility for a light-flavor jet in each event to be mistagged, 
and the probability for a heavy-flavor jet in QCD multijet 
events to pass the CDF 6-tag requirements, respectively. 

Simulated event samples are used to determine de- 
tector acceptance and reconstruction efficiency, estimate 
the contribution from the rest of the SM backgrounds 
with heavy-flavor jets in the final state, and compute 
the number of expected SUSY signal events. Samples 
of simulated it and diboson (M^H^, WZ, and ZZ) pro- 
cesses are generated using the pythia 6.216 [15| Monte 
Carlo generator with Tune A [l6|, and normalized to 
next-to- leading order (NLO) predictions [17|, |l8[. Sam- 
ples of simulated Z/7*(— > ^+^~)-|-jets (Z — e,^, r), Z{-^ 
vv)+]eis, and W{-^ lv)+]ets events with light- and 
heavy-fiavor jets are generated using the ALPGEN 2.1 pro- 
gram [19| interfaced with the parton-shower model from 
PYTHIA. The normalization of the boson plus jets heavy- 
flavor samples includes an additional multiplicative fac- 
tor kHF = 1.4 ± 0.4 which brings the predicted light- 
to heavy-flavor relative contributions in agreement with 
data [ll|. The complete .Z/7*-t-jets and PF-l-jets samples 
are then normalized to the measured Z and W inclu- 
sive cross sections [20|. Finally, samples of single top 
events are produced using the MADEVENT generator 21 1 
and normalized according to NLO predictions 22] ■ The 
SUSY samples are generated in the framework of MSSM 
using PYTHIA. Masses hi and Xi &re flxed, and 6i is only 
allowed to decay via bi — >■ bxi channel. A total of 106 dif- 
ferent samples have been generated with sbottom mass 
M~ in the range between 80 GeV/c^ and 280 GeV/c^ 

and neutralino mass M~a up to 100 GeM /c? . The sam- 

pies are normalized to NLO predictions, as implemented 
in PR0SPIN02 [H], using CTEQ6.6 [H] parton distribu- 
tion functions (PDFs) and renormalization and factor- 
ization scales set to M~ . The production cross section 
depends on M- and decreases from 50 pb to 0.01 pb 
as the sbottom mass increases. The Monte Carlo gen- 
erated events are passed through a full CDF II detector 
simulation (based on GEANt3 [23| and GFLASH |20|) and 
reconstructed and analyzed with the same analysis chain 
as for the data. 

An optimization is performed with the aim to maxi- 
mize the sensitivity to a SUSY signal across the &i — Xi 
mass plane. For each of the 106 signal samples con- 
sidered, the procedure maximizes S/vB, where S de- 
notes the number of SUSY events and B is the total 
SM background. As the difference AM = M~ - M~o 

° bl Xi 



and A(p{]^rp — jet(2)) for the second leading jet, and Ht, 

defined as Ht — J2i ^t (i=l"2)j increase. The results 
from the different SUSY samples are combined to define 
two single sets of thresholds (see Table 1) that maximize 
the search sensitivity in the widest range of sbottom and 
neutralino masses at low AM {AM < 90 GeV/c^) and 
high AM {AM > 90 GeV/c^), respectively. As an ex- 
ample, for M~o — 70 GeV/c^ and M~ in the range from 

Xi bl 

150 GeV/c^ to 250 GeV/c^, values for S/Vb between 10 
and 2.5 are obtained, corresponding to SUSY selection 
efficiencies of 3% to 10%. 



Optimal thresliolds | 


Vt 


Ht + J^t £^t''" ^t'^'" A</.(^j,-jet(2)) 




GeV radians 


low AAf 60 


165 80 25 0.7 


high AM 80 


300 90 40 0.7 



ijct(l) 



increases, the optimal thresholds on ^j,, E; 



ijot(l) 

T ' 



E; 



ijot(2) 



TABLE I: Optimized thresholds on $j,, Ht + $t, E) 
E^^'^'^\ and A(J){$t - jet(2)) for each AM region. 



A number of control samples in data is considered to 
test the validity of the SM background predictions. Sam- 
ples dominated by QCD multijet events or mistags are 
obtained by reversing the requirements on either Iflj, or 
Erp , or using a non-overlapping sample selected with 
looser 6-tagging requirements. Control samples domi- 
nated by Z/7*-t-jets, W-|-jets, and top-quark processes, 
with highly energetic leptons in the final state, are ob- 
tained after requiring either £^xcm/^T -^ ^.9 for at least 
one of the jets, or selecting events with isolated tracks 
with pT.track > 10 GcV/c. Agreement is observed be- 
tween the data and the SM predictions in each of the 
control samples. 

A detailed study of systematic uncertainties was car- 
ried out [10[. The uncertainty on the SM background 
predictions is dominated by the determination of the b- 
jet mistag rates, which propagates into an uncertainty 
in the SM prediction between 13% and 11% as AM in- 
creases. The uncertainty on the kup value applied to 
the boson plus jets heavy-flavor samples translates into 
an 11% uncertainty in the SM predictions. The depen- 
dence on the amount of initial-state and final-state ra- 
diation in the Monte Carlo generated samples for top, 
boson plus jets, and diboson contributions introduces a 
6% uncertainty on the SM predictions. Other sources of 
uncertainty on the predicted SM background are: a 3% 
uncertainty due to the determination of the 5-tagging 
efficiency in the Monte Carlo simulated samples; a 3% 
uncertainty from the uncertainties on the absolute nor- 
malization of the top quark, diboson, W , and Z Monte 
Carlo generated processes; and a 2.5% uncertainty from 
the knowledge of the jet energy scale [9|. In addition, 
uncertainties related to trigger efficiency and the heavy- 
flavor QCD multijet background contribute less than 1% 



to the final uncertainty. The total systematic uncertainty 
on the SM predictions varies between 19% and 18% as 
AM increases. In the case of the MSSM signal, various 
sources of uncertainty on the predicted cross sections at 
NLO, as determined using prospino2, are considered: 
the uncertainty due to PDFs is computed using the Hes- 
sian method [27] and translates into a 12% uncertainty 
on the absolute predictions; variations of the renormal- 
ization and factorization scale by a factor of two change 
the theoretical cross sections by about 26%. Uncertain- 
ties on the amount of initial- and final-state gluon radi- 
ation in the MSSM Monte Carlo generated samples in- 
troduce a 10% uncertainty on the signal yields. The 3% 
uncertainty on the absolute jet energy scale translates 
into a 9% to 14% uncertainty on the MSSM predictions. 
Other sources of uncertainty include: a 4% uncertainty 
due to the determination of the 6-tagging efficiency, and 
a 2% to 1% uncertainty due to the uncertainty on the 
trigger efficiency. The total systematic uncertainty on 
the MSSM signal yields varies between 30% and 32% as 
AM increases. Finally, an additional 6% uncertainty on 
the quoted total integrated luminosity is also taken into 
account in both SM background and SUSY signal pre- 
dictions. 

Figure [1] shows the measured If^rp and Ht + $t dis- 
tributions compared to the SM predictions after all final 
selection criteria are applied. For illustrative purposes, 
the figure indicates the impact of two given MSSM sce- 
narios. The data are in agreement with the SM predic- 
tions within uncertainties for each of the two analyses at 
low and high AM. In Table 2, the observed number of 
events and the SM predictions are presented in each case. 
A global x^ test applied to all data points in Fig. [U and 
including correlations between systematic uncertainties, 
gives a 30% probability for data to be consistent with the 
SM. 



- SM 

' Total Syst. Uncertainty 



- CDF Runll DATA (L=2.65 W) 
■ SM + MSSM 



(2.65 fb-^) 


low AM 


high AM 


mistags 


51.4 ±18.2 


18.5 ±5.5 


QCD jets 


7.6 ±1.9 


1.6 ±0.2 


top 


21.2 ±3.4 


7.8 ±1.3 


Z — > vu+]c\,s 


27.7 ±8.8 


10.9 ±3.5 


Z/'-i* -)■ i+/~-|-jets 


0.5 ±0.2 


0.11 ±0.04 


W — >■ lu+jets 


22.3 ±7.3 


7.3 ±2.4 


WW, WZ, ZZ 


3.1 ±0.5 


1.4 ±0.2 


SM prediction 


133.8 ±26.4 


47.6 ±8.8 


Events in data 


139 


38 



TABLE II: Number of events in data for the two analyses 
compared to SM predictions, including statistical and sys- 
tematic uncertainties summed in quadrature. 

The results are translated into 95% confidence level 
(C.L.) upper limits on the cross section for sbottom pair 
production at given sbottom and neutralino masses, us- 
ing a Bayesian approach [28| and including statistical and 
systematic uncertainties. For the latter, correlations be- 
tween systematic uncertainties on signal efficiencies and 
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FIG. 1: Measured ]^j, and Ht + ^-j^ distributions (black dots) 
for low- AAf (top) and high- AM analyses (bottom) , compared 
to the SM predictions (solid lines) and the SM+MSSM pre- 
dictions (dashed lines). The shaded bands show the total 
systematic uncertainty on the SM predictions. 



background predictions are taken into account. For each 
MSSM point considered, observed and expected limits 
are computed separately for both low- and high-AAf 
analyses, and the one with the best expected limit is 
adopted as the nominal result. Cross sections above 0.1 
pb are excluded at 95% C.L. for the range of sbottom 
masses considered. Similarly, the observed numbers of 
events in data are translated into 95% C.L. upper limits 
for sbottom and neutralino masses, for which the uncer- 
tainties on the theoretical cross sections are also included 
in the limit calculation, and where both analyses are com- 
bined in the same way as for the cross section limits. Fig- 
ure [2] shows the expected and observed exclusion regions 
in the sbottom-neutralino mass plane. For the MSSM 
scenario considered, sbottom masses up to 230 GcV /c? 
are excluded at 95% C.L. for neutralino masses below 
70 GeV/c^. This analysis extends the previous CDF lim- 
its |4| on the sbottom mass by more than 40 GeV/c^. 

In summary, we report results on a search for sbottom 
pair production in pp collisions at ^/s — 1.96 TeV, based 
on 2.65 fb~^ of CDF Run II data. The events are se- 
lected with large $Jrp and two energetic jets in the final 
state, and at least one jet is required to originate from a 
b quark. The measurements are in agreement with SM 
predictions for backgrounds. The results are translated 
into 95% C.L. upper limits on production cross sections 
and sbottom and neutralino masses in a given MSSM sce- 
nario for which the exclusive decay &i — >■ &Xi is assumed, 
and significantly extend previous CDF results. 
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FIG. 2: Exclusion plane at 95 % C.L. as a function of sbottom 
and neutralino masses. The observed and expected upper 
limits from this analysis are compared to previous results from 
CDF and DO experiments at the Tevatron in Run II Q [H, 
and from LEP 29] experiments at CERN with squark mixing 
angle 6 = 0°. The hatched area indicates the kinematically 
prohibited region in the plane. 
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